Nonlinear instability of wormholes supported by exotic dust and a magnetic field 
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Recently, spherically symmetric, static wormholes supported by exotic dust and a radial mag- 
netic field have been derived and argued to be stable with respect to linear radial fluctuations. In 
this report we point out that these wormholes are unstable due to the formation of shell-crossing 
singularities when the nonlinearities of the theory are taken into account. 
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I. INTRODUCTION 

The most popular wormhole model fT] is described by the spacetime {M, g) where 

M = R^xS^, g = -dt^ + dx^ + {x^ + b^){dd^ +sm^^dip^), (1) 

{t,x) denoting Cartesian coordinates on R'^ and i^a) the standard angular coordinates on S^. Here, 6 > is the 
areal radius of the throat which consists of the cylinder x = 0. In order to satisfy Einstein's field equations, this 
solution needs to be supported by "exotic" matter, that is, matter fields which violate the null energy condition. 
Several such matter models have been proposed which admit the solution ([T]). For example, it has been shown 
that this solution may be supported by a ghost scalar field (a massless scalar field whose kinetic energy has a reversed 
sign). Also, the wormhole ([1]) may be_sourced by an anisotropic perfect fluid with an appropriate equation of state 
0. More recently, it has been shown 6] that this solution can also be supported by a combination of exotic dust (a 
pressureless fluid with a negative energy density) and a radial magnetic held. 

Once a suitable matter model which admits the wormhole configuration ([T]) has been found, a relevant question is 
whether or not the solution is stable with respect to fluctuations of the metric and matter fields. Clearly, the stability 
properties depend on the specific matter model. In the case of black holes, for instance, the magnetically charged 
Reissner-Nordstrom black hole is linearly stable in Einstein-Maxwell theory [7|; however, it is unstable in Einstein- 
Yang- Mills-Higgs theory Q if the horizon is sufficiently small. In fact, a similar situation occurs for the Schwarzschild 
black hole when it is embedded in a spacetime with one extra compact dimension to form a black string. While 
the Schwarzschild solution is linearly stable in vacuum general relativity Q the corresponding black string is linearly 
unstable if the circumference of the extra dimension is large enough [lO| . Concerning the wormhole model ([1} , it has 
been shown to be unstable with respect to linear, radial perturbations [ll| when coupled to the ghost scalar field. 
Numerical simulations in the nonlinear, spherically symmetric case indicate [T2l - [l^ that the wormhole either rapidly 
expands or collapses to a Schwarzschild black hole. Attempts have been made to stabilize the wormhole by adding an 
electric or magnetic charge Q^, resulting in charged generalizations of ([1]). Although the addition of charge reduces 
the time scale of the instability, linearly stable wormholes have not been found in the scalar field model. 

In a recent article flGj the linear stability of the wormhole ([T]) has been analyzed in the matter model put forward 
in with the conclusion of stability. In this report, we point out that the wormhole ([IJ is unstable when nonlinear, 
radial fluctuations of the metric and the exotic dust are taken into account. The reason for the instability is that 
for a large class of initial data shell-crossing singularities form and therefore curvature invariants of the perturbed 
spacetime blow up in finite time. 



II. THE MODEL 

The model considered in [(| consists of a spherically symmetric, self-gravitating dust cloud coupled to a radial 
magnetic field. In terms of a Lagrangian coordinate x which parametrizes the dust shells and proper time r measured 
by an observer comoving with such a shell, the key equation is a master equation^ for the areal radius r(r, x) which. 



^ See the Appendix for a derivation. 



2 



for each fixed shell x, has the form of a one-dimensional mechanical system, 

1 



-f{T,xy + V{r{T,x),x) ^ E{x), 



with potential 



y(r,x) 



H{x) 



2r2 



(2) 



(3) 



Here, the energy function E{x) is determined by the initial data ro(a;) := r(0,x) and vo{x) := r(0, x) for r{T,x) at 
r = 0, say, Q is the magnetic charge and /Lt(a;) is the mass function for the dust, which, in terms of the initial density 
po{x) is 



H{x) ^ fioo- 47rG / Po{y)ro{yfr'Q{y)dy, 



(4) 



fioo representing the mass at a: — >■ oo and G Newton's constant. Once the function r(T, x) has been determined, the 
metric g, the electromagnetic field tensor F, the four-velocity u and the density p are given by 



g = -dT^ 



r'{T,xf 
1 + 2E{x 
F = Q d-d A s'm-d dip, 
d 



dx^ +r{T, xf(di)^+sm^'ddLp^), 



p{t,x) = po{x) 



r{T, x) J r' [t, x) 



(5) 
(6) 
(7) 



Here and in the following, a dot and a prime denote partial differentiation with respect to r and x, respectively. 
Notice that there is a freedom in reparametrizing the dust shells, x ^ fix). In the Tolman-Bondi models (see, for 
instance, Ref. 17]) this freedom is usually fixed by labeling each shell by its initial areal radius, ro{x) = x. Since here 
we are interested in describing wormhole spacetimes with one throat, we set ro(x) = y/x'^ + b'^ instead, where 6 > 
is the throat's initial areal radius. 

Static solutions are obtained if each dust shell x is in equilibrium with respect to the potential V{r,x). This is 
possible if both and p{x) are strictly positive, in which case there is a potential well with global minimum at 
r = Q'^/p{x) corresponding to the energy E{x) = —Q^/{2r^). This, together with r{x) = ro{x) = ^x'^ +b'^ leads to 
the static wormhole solution Q with 6 = \Q\, negative density 



Po{x) = - 



(8) 



and zero asymptotic mass pac — 0. 



III. PERTURBATION ANALYSIS 



Now let us consider a nonlinear perturbation of the static solution described in the previous section for a fixed value 
of the magnetic charge Q. Its Cauchy evolution is determined by the initial data for the functions ro{x), vq{x) and 
Po (x) and the value for the asymptotic mass p^o ■ Keeping the labeling of the dust shells such that vq (x) = \J x^ -t- 6^, 
this amounts in perturbing vq{x), po{x) and poc from their equilibrium values vq = 0, po given by Eq. ([S]) and poc = 0, 
respectively. The perturbations of po and pa^ give rise to a new mass function p{x) according to Eq. For the 
following, we consider a rather large class of nonlinear perturbations which are "small" in the following sense: 

(i) The perturbed mass function p{x) is strictly positive for all x G M. 

(ii) The total energy E{x) is strictly negative for all a; G M. 

Condition (i) implies that V{-,x) describes a potential well with negative global minimum at rmin(x) = Q'^/p{x). 
Condition (ii) means that the trajectories of the system ([2]) are bounded; hence the areal radius of each dust shell 
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undergoes periodic oscillations about rmin{x)- If < ri{x) < r2{x) denote the zeroes of the function r n- E{x) — V{r, x) 
for each fixed x, the corresponding period is 

r2 (2;) 

T(x) = 2 [ f9) 

ri{x) 

Since the static wormhole solution satisfies the conditions (i) and (ii) it is not difficult to construct initial data with 
these properties. For example, one could leave po and fico unperturbed which also leaves the potential unperturbed 
and choose \vo\ small enough such that Vo{x)'^ < —2V{rQ{x),x) = Q^/ro(x)^ for all x S R. 

Since small perturbations fulfilling the conditions (i) and (ii) behave like a one-dimensional particle in a potential 
well one might expect the static wormhole ([1]) to be stable under perturbations. However, a problem arises when 
different shells of dust touch or cross each other, in which case the density function p(t, x) diverges. This occurs at 
points where the function vir, x) := r'{T, a;)/rg(x) vanishes, see Eq. Q. The time evolution of vir, x) is governed by 
the following equation which may be obtained by differentiating Eq. ^ twice, 

i)[T,x)+u(T,xfy{T,x)^-^T,Gpo{x)(^^\ , (10) 

\r{T, x)J 



where 



^[r,x)^—^{r{T,x),x) = -' 
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dr^ ' ' r(T, x)^ r{T, x) 

Since ^"(0,2;) = 1, i'(0,a:) — Wg(a;)/rQ(x), shell-crossing is absent for small enough time provided that the initial 
function 

i){0,x) = ^(x) 
ar 

is regular. In particular, this requires the initial velocity profile vq{x) to have a critical point at the throat. Otherwise, 
the time evolution of the wormhole is not defined since in this case 



p{0,x) = -po{x) 



2—- -I- u[0,x) 
rQ[x) 



diverges. The regularity of i>(0, x) also implies the regularity of the radial part of the metric. 



whose Gaussian curvature is u~^v, for small enough time. However, there is a coordinate singularity at x = unless 
^o(^)^/ (1 + 2E{x)) is finite and different from zero at a: = 0. This yields an additional restriction on the initial energy 
E{x). In order to shed light on this condition, we assume for simplicity that po and poo are unperturbed and set 

vo{x) = a{x)\Q\/rQ{x) with some smooth fimction a{x) satisfying a(x)'^ < 1 for all a; G M. Then, 

l + 2E{x)^ ^—:r{x^ + Q^a{xf). (12) 

Since r'Q{x) — x/vq^x), the coordinate singularity is avoided if a{x) is proportional to x for small \x\. Together with 
the regularity of i'{0,x) this means that the function a{x) should satisfy 

a(0) = 0, q'(0) = 0, a{x)^ < 1, xeR. (13) 

In the following, we assume that the above conditions on a{x) are satisfied in order for the coordinate x to be 
globally defined on the real axis and in order to avoid shell-crossing singularities for arbitrary small times. To analyze 
the occurrence of shell-crossing singularities at later times, we differentiate the equation 

r{T + nT{x),x) — r{T,x), n = 0, 1, 2, 3, r, xgM, 

and obtain 

T'(x) 

v{t + nT(x),x) — v{t, x) — nr(T, x) , , . . (14) 

raix) 
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The evaluation of this equation at a turning point t = ti{x), where r{Ti{x),x) = 0, yields 

i'{ti{x) + nT{x), x) — iy{Ti{x), x)^ n = 0, 1, 2, 3, x G M, 

while for r = we obtain 



v{nT{x),x) 



1 



, ,T'{x) 
n'Vo[x) 

r'aix) 



0,1,2,3, 



X £ 



(15) 



(16) 



Therefore, for each fixed shell x, the field v{t,x) oscillates between the two values given in the right-hand side of 
Eqs. ([T5|) and ([T6l) . respectively. If vq{x)T' {x) / r'f^{x) ^ 0, this signifies that the oscillation's amplitude grows without 
bound for t oo. Furthermore, since v{Q, x) ~ 1, Eq. ([T6l) also implies that for vq{x)T' {x) / r'f^{x) > a shell-crossing 
singularity must occur in finite proper time. It is a simple matter to construct initial data satisfying this condition 
everywhere away from the throat. For example, leaving po and /ioo unperturbed and setting vo{x) = a(x)|Q|/ro(x) 
with the function a{x) satisfying the restrictions ()13p gives 



= 4^a(a;)(l - a{xf)-'''^ 



1 -a(x)2 + -ro{x)' 



a{x)a'{x) 



If ao{x) is any smooth function which satisfies the conditions (IT3|) . is strictly positive for x 0, and constant outside 
a compact interval, then the family of functions a a (a;) = Xao{x) with A > small enough is easily seen to satisfy 
vo{x)T'{x)/r'fy{x) > for all x ^ 0. We conclude that a large class of nonlinear perturbations leads to the formation 
of shell-crossing singularities. As a consequence, the wormhole is unstable. 

We shall finish this report by making a few comments regarding the behavior of linear perturbations on the static 
wormhole background defined in the previous section. If we set r{T,x) = ro(x) -I- dr{T,x), differentiate Eq. ([2]) with 
respect to r, we obtain, to linear order in (5r(r, x), 

d'^Sr(T,x) / x?c- / X „ / ^ IQI 



9t2 



rlix)'- 



and thus. 



Sr{T, x) ~ Svo{x) 



sin(a;o(a;)T) 



(18) 



ujo{x) 

where we have assumed that rg, po and /^oo are unperturbed. For the following, we also assume that the linearized 
initial velocity 6va(x) is compactly supported away from the throat at a; = 0. The induced linear density perturbation 
Sp{x, t) is 



Sp{T,x) = pq{x) 



2Sr{T, x) Sr (r, x) 



(19) 



ro(x) r^ix) 

Since each shell oscillates with its own frequency iL)o{x), it is seen that 5p{T, x) exhibits linear growth in time (see also 
the discussion in fl6l|). In terms of the gauge-invariant combination 



this yields 



J(t,x) 



J(t, x) :— 5p{T, x) 



ro{xYpo{x) ( dr{T,x) 



ro(x)2 



-Sr{T,x), 
r'oix) 

ro{x)^Po{x) ( si n(a;o(x)T) 5vo[x) 
ujq{x) ro(x)2 



(20) 



Therefore, the linear growth in t cannot be transformed away by an infinitesimal coordinate transformation. As we 
have shown above, in the nonlinear case, this growth leads to blowup in finite proper time. 



IV. CONCLUSIONS 



The result of this report is that the wormhole ((T)), when supported by exotic dust and a radial magnetic field, is 
unstable due to the formation of shell-crossing singularities. Although shell-crossing singularities have been argued 
to be "mild" in the sense that it is often possible to construct a C^-extension of the metric [3], nevertheless the 
density function and the Ricci scalar blow up. Even if one is willing to accept such singularities it is worth pointing 
out the following: according to the arguments presented below Eq. (|16p . the extended spacetime possesses the bizarre 
property that for large t, the norm of the vector field d^, measuring the normal geodesic deviation between neighboring 
dust shells, undergoes wild oscillations. 
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Appendix: Derivation of the master equation 



The radial components of the Einstein equations for the metric 

g = ^dr' + dx^ + r{T, xf[d§^ + sin^ § V), 

the four- velocity u—-^ and the electromagnetic field Y ^ Q d-d f\ sin d dip read 

2 r 2m' 

G% = 7f'^-^--^-8-Gp, (21) 
r- 2.. 2m Q2 

G% = -^r' = 0, (23) 
r 1 

where G"/3 are the components of the Einstein tensor and TO(r, R) is the Misner- Sharp mass function [l9j which is 
defined by 

l- — = g{dr, dr) = -r2 ^ . (24) 
r 

These equations are solved as follows: first, Eq. ([25)1 immediately implies that F = F(x) is a function of x only. 
Next, differentiating Eq. ([M)) with respect to r and comparing with Eq. yields m = (5^r/(2r^) which may be 
integrated to obtain 

m(T,x) =Ai(x)-gV(2r(T,x)), (25) 

with /i(a;) representing the dust mass function. Finally, using this in Eq. (PT|) gives fj! — AirGpr^r' . As a consequence, 
pr^r' must be independent of t which implies Eq. ([7|), and Eq. Q is obtained after integrating in x. The key equation 
([2]) that governs the dynamics of the system is obtained from Eq. (|24)) after setting F(a;)^ — 1 -|-2£'(a;) and taking into 
account the relation (1251) between the mass functions m and /i. The remaining components of the Einstein equations 
and the mass conservation law are automatically satisfied as a consequence of the previous equations and Bianchi's 
identities. 
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